Some Mathematical Models for ELM Signal 
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There is no wide accepted theory for ELM (Edge Localized Mode) yet. Some fusion people feel 
that we may never get a final theory for ELM and H-mode, since which are too complicated (also 
related to the unsolved turbulence problem) and with at least three time scales. The only way out 
is using models. (This is analogous to that we believe quantum mechanics can explain chemistry 
and biology, but no one can calculate DNA structure from Schrodinger equation directly.) This 
manuscript gives some possible mathematical approaches to it. I should declare that these are just 
math toys for me yet. They may inspire to good understandings of ELM and H-mode, may not. 
Useful or useless, I don't know. One need not take too much care of it. Just for fun and enjoying 
different interesting ideas. 



I. INTRODUCTION 



II. ORDINARY DIFFERENTIAL EQUATION 
(ODE) 



Typical ELM signals are shown in FIGQ] Here, we try 
to model type-I and type-Ill ELM signals of D a . 



Pulse No: 52994 



15 
141 



5- 




1.6 

IJ 



0.4 



1.5 
ri. 1.0 
cf 0.5 



T. (KfiV) 




L-H Transition 

i Type-I H 
L-Mode ' ELMs 



rm — rrm" 
Type-I ELMs 



ELM-free 
,H-Mode 



17.0 



17.5 



-t~ 




1B.0 ^8.5 
Time (s) 



19.0 



FIG. 1: JET typical ELM signals ([Perez2004, thesis]). 



We just focus on the signal itself and ignore all the 
other information behind, i.e., we just discuss how to use 
(as simple as possible) equations or other approaches to 
reproduce the shape of these signals. It is qualitative or 
at most semi-quantitative. 
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This is a wide used approach for modeling. For ex- 
ample, the famous prey-predator model for fishbone and 
drift wave-zonal flow system. One may also extend it 
to PDE to contain the spatial information. Indeed, this 
way has been used to model ELM and L-H transit in 
literatures, e.g., [Diamondl994], [Itoh.1991, 1993, 1999]. 

If we assume that the dynamic of signal is not explicit 
depend on time t, then one equation y = f(y) is not 
enough because that y is not single valued with y. So, 
at least, we must use second order equation or use two 
first order equations. While, one second order equation 
is just a special case of two first order equations. 



A. Example 1 

In this example, we use the second option, i.e., two first 
order equations, which is also used in [Diamondl994]. 
The equations are 



m = a m ~ hyl - b 2 yiV2, 
V2 = -a 2 y2 + b 3 yiy 2 . 

A result is shown in FIG. [2j 
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ODE modeling for ELM (1), a 1= 1, a 2 =0.5, b 1= 0, b 2 =10, b g =1 
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FIG. 2: ODE for modeling ELM, example 1. 

The physical explanations of the equations and param- 
eters can be found in [Diamondl994], however, where 
2/1 = E,y 2 = U. 
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B. Example 2 

Using one second order equation, e.g., 

y = a(l - by 2 )(y +c)-y-ac, 



DDE for modeling ELM, example 1 



(2) 



which is modified from stiff system equation, a result is 
shown in FIG.rj 

The shape is very similar to the PDE result in 
[Itohl993]. But, an apparent drawback of ([2} is that 
the signal y is not always positive here. 



ODE for modeling ELM (2), a=5, b=0.4, c=0.12 
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FIG. 3: ODE for modeling ELM, example 2. 



III. DELAY DIFFERENTIAL EQUATION (DDE) 

DDEs contain derivatives which depend on previous 
time. We guess some equations here. The equations are 
from or modified from [Shampine2000] . There is also a 
type of DDE has sawtooth solutions (see e.g., [Mallct- 
Paret2011]), which may be also suitable for model the 
unsolved sawteeth phenomena in magnetic confined fu- 
sion study. 



A. Example 1 



The equation is 



y'(t) = -Xy(t-l)(l + y(t)), 
y(t)=t,t<0. 

Results are shown in FIG. [4] 



(3) 



B. Example 2 



Equation is 



y (t- 0.74) 



y(t) = 19,t< 0;y(0) = 19.001. 



y'(t) = ry(t) 1 



(4) 



which is original to model four-year cycle of the popula- 
tion of lemmings. 

The result is shown in FIG. [3] 
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FIG. 4: DDE for modeling ELM, example 1. 



DDE for modeling ELM, example 2 



Phase Plane, r=3.5, m=19 




IV. 



FIG. 5: DDE for modeling ELM, example 2. 



DIFFERENTIAL EQUATION + MONTE 
CARLO (DEMC) 



This is a more intuitive way to model ELM, which com- 
bines different time scales using a more acceptable ap- 
proach. The burst of ELM signal is suggested related to 
peeling-ballooning mode in standard ELM model ([Con- 
norl998]). So, we give two time scale: transport time 
scale and ballooning mode time scale. 

To simplify the transport process, the transport time 
scale is modeled using 



> - 

at tt 



(5) 



When pressure gradient exceeds marginal value (using 
f c i to represent it), the ballooning mode occurs. The 
signals will grow exponentially 



= 7bm/, 



(6) 



with the characteristic time tbm = 1/7-BM- 

When f{t) exceeds a marginal value f C 2, the plasma 
crash. And then, go to next cycle. We also add some ran- 
domicities to / cl and / c2 to make the whole signal more 
like the actual experimental signal. A result is shown in 
FIG.E1 
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DE+MC for modeling ELM, ^=1 .062, x 2 =0.001 
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FIG. 6: DEMC modeling for ELM. 

This simple DEMC model captures some fundamental 
properties of integrated simulations, such as [Chang2008] 
and [Lonnroth2009]. 



V. SUMMARY AND MORE 



Which one is more reasonable? At present, I choose 
DEMC > ODE > DDE. While, maybe none of them can 
be final explanation of ELM or H-mode. 

There is another math tool called Fractional Differen- 
tial Equation (FDE) , which has been used for anomalous 
transport/diffusion problems in some areas. For plasma 
community, FDE is not wide accepted yet, but it is pos- 
sible in the future. FDE is hard to give periodic solutions 
as ELM at present. So, I haven't given examples. 
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